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Abstract 

This article addresses the important problem of electrodynamic coupling of 
3<1 stellar coronal loops to underlying velocity fields. A rigorous analysis 
has revealed that the physics can be represented by a simple yet QXivJ,V<tZZYlt LRC 
circuit analogue. This analogue points to the existence of global structure 
oscillations which resonantly excite Internal field line oscillations at a spatial 
resonance within the coronal loop. Although the width of this spatial resonance 
as well as the induced currents and coronal velocity f ield^ explicitly depend upon 
viscosity and resistivity, the ACAoncint form of the generalized electrodynafi'i.c 
heating function is virtually >Lndzpmdmt of irreversibilities. This is a classic 
feature of high quality resonators that are externally driven by a broad-band 
source of spectral power. 

The major results of this article are: 

(1) The heating function, E^j, field-aligned electron current, and cross-field ion 

polarization current are explicit functions of the 3^1 velocity field's spectral 

power Junction, e.g. , E„=6.28xlO^T^^f/il2^^) (^loop^ 

The essential feature of a resonant heating mechanism is that magnetic loops 
with different lengths, and hence different global resonance frequencies, 
are heated at a rate that critically depends upon the amount of spectral 

il 
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power found at the resonance frequt>ncy, V 

(2) This heating function results In the following IjQpliclt equations for the 

naxlmim teoperature, T, and base pressure, P: 

4 r , photo 

T - 3.36 X 10 exp (0.36f^^^/fp) < l/2pv‘ 

1 3g X 10^ / \r -i photo! 

F - j Ktf ^ \ J 

which are In remarkable agreement with observations of solar coronal loops. | 

(3) Solar applications require nonlinear modifications that do not alter (1) and j 

(2) above but which dramatically increase the cross-field size of the | 

j 

coronal dissipation site from kilometers to thoUACLndi of kilometers. Preferential 
viscous heating of the ions is also predicted as well as a 10-16 km/sec coronal 
velocity field. 

j 

i 

i 

1 

] 

j 



ill 


RESONANT ELECTRODYNAMIC HEATING OF STELIAR CORONAL LOOPS: 


AN LRC CIRCUIT ANALOGUE 
I . Introduction 

One of the most outitanding probloma faced by solar and stellar physicists 
is an understanding of the cause and effect relationship between inner (e.g», 
photospheric) and outer (e.g. , coronal) atmospheric dynamics. Central to this 
problem of global coupling within a stellar atmosphere is the role playtd by 
magnetic fields. Magnetic fields typically thread the inner and outer atmospheree 
of a star and thus one cannot disregard the potentially important role that 
associated electrodynamic processes could play in effecting an "electrodynamic 
coupling" of these two regions. 

A useful parameter that identifies the dominance of "electrodynamic coupling" 

over, for example, mechanical coupling of the inner and outer atmosphere by 

2 2 

acoustic waves, is the plasma beta, p = v /v (v « sound speed; v Alfven 

S A 8 A 

speed). Specifically, electrodynamic coupling within a stellar atmosphere 

dominates mechanical coupling whenever the plasma beta is less than unity within 

the outer atmosphere, l.e., p < 1, and of the order of unity within the 

inner atmosphere, i.e., ^ Thus, mechanical dynamics such as convective 

and/or differential-rotation velocity fields within the ^ T atmosphere 

can couple to and drive casually related phenomena, such as heating, within the 

P <1 outer atmosphere through various electrodynamic processes associated 

outer 

with the Interconnecting magnetic field. This is in contrast to a mechanically 
coupled stellar atmosphere which would require a p 1 outer atmosphere. 


? 

It now widely believed that the 0 > 1 inner (photosphere) end 3 < 1 
outer (corona) atmosphere of the Sun are electrodynamically coupled since 
extrapolations of the observed photospheric magnetic energy into the corona 
exceeds the in*situ thermodynamic energy (i«e., S < 1). Furthermorcy it 
appears that the solar corona is typically structured, comprising a variety 
of closed, loop-like regions of enhanced radiative output, i.e., coronal loops 
(c.f., Uithbroe and Noyes, 1977; Vaiana and Rosner, 1978). In this article 
I will focus upon the electrodynamic coupling of solar coronal loops to their 
underlying driver -- the S ^ 1 photospheric velocity field. However, this by 
no means restricts the applicability of the following analysis t<a the solar 
setting. Indeed, the results of this article are generally applicable to any 
stellar atmosphere and for that matter any cosmic plasma that satisfies the 
conditions on the plasma 8 as described above. 

The specific goals of this article are to: 

(i) Determine both the resistive response (i.e., an electrodynamic heating function, 
and the reactive response (i.e., the Induced currents and velocity fields) of 
magnetic loops in terms of the spectral power of an underlying 3^1 velocity field 
such as the solar photosphere. 

(li) Utilize the heating functi 'n, E„, in an energy balance model whicli yields the 

H 

resulting thermal state (i.e., maximum temperature and base pressure) of a loop* 
also in terms of the spectral power of the 6^1 velocity field. 
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With regard to (i), there have bean a number of attampta at developing a coronal 
heating function, E„ (c,f., reviews by Wentsel, 1978; Hollweg,1979 ; Kuperus, 
lonson and Spicer, 1981). For example the solar physics cosuMmity has proposed 
heating by Alfven waves, fast waves, slow waves, Alfvenic surface waves, 
electrical currents, magnetij^ raconnection •*> a seemingly endless list of 
possibilities. Unfortunately, there is a very serious problem which appears 
to be rapidly spreading throughout the solar and into the stellar comnunities . 
this problem is associated with the fact that each of the above mentioned 
mechanisms are being thought of as unique in their own right. This, however, 
is simply not true since they are all based upon a unifying foundation of 
electrodynamics whose principle ingredients are capacitance -- the ability of 
a magnetoplasma to store electric and kinetic energy. Inductance -- the ability 
to store magnetic energy, and resistance -- the ability to convert the above 
electrodynamic energy into thermodynamic end products such as heat. Another more 
practical problem with the plethora of proposed heating mechanisms is that they 
do not result in a heaulng function which is explicitly correlated with the 
properties of the 8^1 driver of the corona, viz., the underlying photospherlc 
velocity field. Thus, it has been virtually impossible to develop a complete 
thermodynamic model that allows us to appreciate coronal loop heating in terms 
of synchronously observed mechanical activity within the photosphere. 

Although some authors have invoked electrodynamic circuit analogues In an 
attempt to consolidate the many different aspects of electrodynamics into a 
simple formalism, their approach has for the most part been somewhat too 
phenomenological as well as limited in scope. Specifically, they have primarily 
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addresied flare prcblema In terms of transient discharges (c.f., Alfven, 

1977 and raferences therein). They have not, however, considered continuously 
driven systesw and have also virtually disregarded the critically important capacitive 
properties of the magnetoplasmay Thus, although their basic phenomenological treatment 

deserves coomendation, the concept of electrodynamlc circuit analogues should 
be rationalized by a more complete and self-consistent treatment. Therefore,* 
upon Investigating the Important global cause and effect relationship between a 
coronal loop and Its 0^1 driver, one should keep In mind the possibility that 
the electrodynamics of such a system can be represented by a simple yet O.qiU.vaJLe.i'lt 
driven "LRC* circuit. If this is the case, then a powerful approach 

to the complex problem of "electrodynamlc coupling" In astrophysical plasmal would 
be at hand since the electrodynamics of simple LRC circuits Is well understood. 

Of course, the key Is in the determination of the equivalent L, R, and C values 
which, as we shall see. Is relatively simple to do. 

II. The Physics of Global Electrodynamlc Coupling 

Figure 1 Illustrates both the prototype physical system that will be 
investigated and its equivalent electrodynamlc circuit which will be derived 
below. The physical system comprises both a E < 1 magnetic loop and an 
underlying region of E > 1 velocity fields which electrodynamically couple to 
and drive electrodynamlc activity within the 0 < 1 loop via the Interconnecting 
magnetic field B . The magnetic loop contains plasma at both coronal and 
chromospheric temperatures, the transition occurring along the externally 
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gtncnted potential nagnetic field, The field-aligned scale length of 

the coronal portion of the nuignetic loop will be denoted as f , iiliile the 

cor 

scale length of the chrosMspheric portion will be denoted as 
addition, the field -aligned scale length of the P > 1, photospheric velocity 
fields tdiich interact with the magnetic loop will be denoted by 
cross-field scale size of system will be denoted as which represents the 
diameter of the magnetic loop in the corona, chromosphere and photosphere. 

For simplicity will be assumed to be the same In all three regions -- an 
assumption that will not adversely affect the analysis presented in this article. 

The main field is generated by a primary dynamo that is external to 

the loop's local mechanical driver. This is probably the case for many stellar 

coronal loops in which large scale velocity whose characteristic time scale 

exceeds the loop's Alfven transit time, are responsible for the overall potential 

magnetic field structure. It is the relatively smaller scale 3 > 1 velocity 

fields^ whose characteristic time scale is of the order of the loop's Alfven 

transit time, which drive the loop into a non-equilibrium electrodynamic state 

resulting in the flow of electrical currents along <i.e., force-free) and 

across (i.e«, non-forcefree) the ambient field The equivalent LRC circuit 

is illustrated in hindsight the derivation to follow and models the electro- 

dynamic behavior of the physical system. Specifically, the g ^ 1 velocity fields 

act as a secondary voltage generator, supplying a time dependent emf , which 

drives an equivalent current, I, though an equivalent loop inductance, L, 

resistance, R, , and capacitance, C, as well as through an equivalent photo- 
loop 


spheric resistance, R 


photo* 
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U) VzfUvcUion oi tht Global tltctAX>dynomlcA tqaa^^'i'n 

Magnteohydcodynaiaic syfttms arc inherently nnnllnear. Hoveverp experience 
in handling auch complex eystems has shown that critically important characteristics 
become apparent by appropriate linearization. In thle regard , the following 
equations will be linearized with nonlinear nodifk.wationa being considered later in 
the article: 

4TT hi 

__._7x(vxy (1) 

c ot 


l-l-O ( 2 ) 


% E 

E + Ti j (3) 

c ^ 


A, ** E 

p 3? + PX • S:- * 1,+ P|=+ “7^ 

where 


t = P L - I 

(5) 

T * U • V V 

fSi S ~ 
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p,f - (1 + 
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And where the elertron collieion frequency* end Ion colliilon frequency* 
v^* ere given by; 


V * y V 4* y 

e e-1 e-neutrel enooeloue 


1/2 

''i * <“e^“l^ ''e ^e 


( 10 ) 

(U) 


with the electron-ion collieion frequency* v * given by: 


V - 50,0 (n /T eec’^* 
0-t e 6 


the electron -neutral collision frequency* v given by; 

c-n 


~ -1 
V ■ 1,95 X 10 n T sec 
e-n n c 


( 12 ) 


( 13 ) 


and with the anomalous collision frequency, '’gjjomloyg* electrons with a 

non-thermal level of low frequency (40 < Ui .) electrostatic plasma waves being 

pi 


given by: 


^anomalous 


0 ) ( " ■■■ \ 
\"eVe/ 


( 14 ) 


where W is the energy density of the plasma microturbulence. Such plasma 
microturbulence could develop if the electrical currents exceed various 
threshold conditions (c.f.* Papadopoulos * 1977; Ions on et al,. 1979; Rosner 
et al. , 1978; Hinata, 1980), 

The electrodynamics equation (1) is given in terms of the local current 
density, and electric field, E. In analyzing this equation I will adopt 
a cylindrical coordinate system in which X|| represents a direction parallel to 
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th« Mblant fitld, and in ii«hlch raprasants a diractlon parpandieular to 
tha fiald. Tha parpandieular unit vactor, contains two basaS| viz., Xg 
«dilch corraaponda to the azinuthal diraction and x^ which corraapondi to tha 
radial diraction* In addition* variations in plasoa paraoMtars such as tha 
Alfvan spaad will occur only in tha xy and x^ diractions. 

In order that tha Poynting flux ba directed primarily along tha ambient 
field* J^^* it follows that (7 x E)|| 0 (i*a.* approximately shear disturbances). 

Note that the field aligned perturbation in the magnetic field is not identically 
zero thereby resulting in a small component of the Poynting flux in the x^ 
direction. This is Important since it allows the "interior" of a coronal loop 
to receive electrodynamic energy "exterior" regions of a loop. Note also that 
the condition <7 x E)|| 0 is equivalent to disturbances that are approximately 

incompressible (i.e.* ~ 0) thereby rationalizing the use of an incompressible 

momentum equation (4). Although nonlinearities could convert shear into 
compressible disturbances therefore requiring use of a compressible momentum 
equation (Hollweg, 1981)* it is important to first consider the linearized, 
shear system. In section IX. (e), I will briefly discuss some simple nonlinear 
modifications of the linear analysis. 

» n 

Under the sonstralnt that (7 x E)ii ~ 0 it follows that 7 (7 * E ) - V 

and thus equation (1) can be written in the following component form: 




ot 




E ) - 7 



+ 


0 


(15) 


9 



(16) 


NoClng in hindsight that d/StK » due to a local resonance absorption 

r y 

process that will be discussed later, it follows that E. « E ; and since 

b r 

1 E the J cuvren-: is carried primarily by j . Using Ohm's law (i.e., 
^X ^X ^ 

equation (3)) the electrodynamics is therefore described by, 



Since 7 • j « 0 (i.e., Bj /Bjc Bj,i/Bx„), equations (17) and (18) are redundant, 

* r r !1 II 

I’his redundcncy is primarily a consequence of the (7 x ^ ^ 0 constraint. Thus 

it is necessary to consider only one of the above equations since they both contain 
a complete electrodynamic description. In this regard I will focus upon equation 
(18) , rewritten in the following form. 


2 




BX||Bx Bt 


H(B-l) 


3 

B (VqB /c) 


BxiiBx Bt 

1! r 


(19) 


where H(l-8) is the Heavyside step function and with P as the plasma beta 
(i.e., H(1«P) * 0 for 3 > 1 and ■ 1 for 3 < 1) and where the second term on the 
left hand side of equation (18) has been neglected since it is small compared to the 


other terms. 
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The third term on the left hand tide of equation (19) contributes only 
in P < 1 regions and represents the hydrodynamic reactance of the plasma to 
externally applied emfs associated with the P ^1 velocity fields (l.e., the 
right hand side of equation (19)). The P < 1 reactance terut can readily be 
determined via the force balance equation (4) valid In P < 1 regions, l.e., 


9v 

_ 

P “ * 


•9 . _2 

+ M 7 V . 
c ~x *^x 


( 20 ) 


Using the correspondence d/dt lU) and noting In hindsight that the reactance 
Is maximum when 3/9t <* follows that 




SJl 


4TTv. 


Pr tl 
roapt 


Sx ^5t 

r 


&x| 


( 21 ) 


where the magnetic Prandtl number Fr Is defined as the ratio of momentum 

mag 

diffusion to magnetic diffusion, l.e.. 


Pr 


mag 


4% 

L. 

2 

pc r\ 


, vl/2 / 2 1-1 

(^) 6*t) ■ 


( 22 ) 


Using equation (21)^ the local electrodynamics equation (19) reduces to, 


4TT 






- n (1 + H(1 - P)Pr „) 


/^Tiv \ aj“ii 


^ (VpB^/c) 
5t 


(23) 
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•nd Cite begins to tee the emergence o£ e simple "LRC" circuit enalogue -- the 
first term representing an Inductive reactance, the second term representing 
a Joule and viscous resistance, the third term representing a capactlve reactance 
and the term on the right hand side representing the external source of emf which 
drives the currents. 

Equation (23) describes the local electrodynamics of the syittem. A 
description of the global electrodynamics Is obtained by Integrating equation 
(23) over the volume of the system. In performing this integration it is 
important to note the possible occurrence of a local resonance absorption 
process In the ^ ^ 1 magnetic loop. This process has been discussed In some 
detail by a number of authors (e.g., Hasegawa and Chen, 1976; Kappraff and 
Tataronls, 1977; lonson, 1977, 1978, 1981; Ott e^ al. , 1979). Although most 
of these discussions are strictly applicable for specific plasma parameters 
(e.g., collisionless or collisional analyses), lonson (1981) has , presented 
the concept in a relatively unified manner which I will briefly review for the 
case of a collisional plasma (i.e., the electron pressure E^*«field which supports 
kinetic Alfven waves is small compared to the resistive E-field). 

Noting that there is a cross -field variation in the Alfven speed, dv^/dx^, 

it follows that for a particular driving frequency, ^absorption of electro- 

dynamic energy could occur within a small spatial bandwidth, Ax^, centered about 

a spatial resonance located at - \es* ^ specific Xjj, a spatial 

resonance exists only if v (x * x )d/&x<|iw at some point in the x -direction— 

A T r€ts X 

a condition that is not necessarily satisfied. However, if this resonance 
condition is satisfied, the spatial bandwidth, Ax^, can readily be estimated from 
equation (23) by expanding the Alfven speed about the spatial resonance and 
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noting that the Inductive and capacitive ri.7actancee cancel one another at the 

2 2 2 2 2 2 

spatial resonance, ■ \es* Ijj/c 5t“ * ^^\^*res^^ ^ll^^l! * 

A consequence of matching the inne,'r ' esii^sive solution to the outer reactive 
solution (i.e*,, the ideal ohd solution) Is that the resistive term in equation 
(23) equals the fir^vt order contribution of the expanded capacitive term, i.e.. 


H(l-P)) U ^ 

SSS 

(Ax^)^ 


ATT u? (Ax ) 

r 


2 

ac 


(24) 



(25) 


and uliere the correspondence iTT/(Ax^) has been used. Solving equation 

(24) for Ax yields, 
r 




(1 + Pr H(1 
mag 



Re 

mag 


-1/3 


( 2 ^) 


vhere the magnetic Reynold's number, defined as, 


Re 


mag 


local . 2 
ATT v, £. 

A x_ 

nc £)| 


(27) 


where 3/2Xjj-> iTT/Xjj has been used with Xjj representing the field-aligned extent 
of the local region. It should be noted that Ax^ as defined by equation (26) 
depends upon local plasma parameters such as the resistivity, ri, Alfven speed, 
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V. , the magnetic Pranitl nuoiber* Pr and the £leld>aligned extent of the 
A Mg 

local region, jtj|. It ia alao important to note that equation (26) Implicitly 

* 

aaaumea the exiatence of a apatial reaonance. 

Keeping in mind that can be interpreted aa the local (i.e., along the 
magnetic field), cheracteriatic croaa^field ahear length of the induced current 
and velocity fields, equation (23) can be aimply Integrated over the volume of 
the system. This results in a description of the system's global electrodynamics 
which is represented by the following global electrodynamics equation: 


L R + 

2 global dt C dt 

at 




( 28 ) 


where the equivalent Inductance, L, capacitance, C, resistance, R 
current, I, and driving emf, ^(t) are given by, 


global’ 


H, 


L * 


Tic 


2^ 


(29) 


je :c 
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(30) 


R - , - * R- ^ j. 

global loop photo 


(31) 




^ 11 (1 + Pr ) X||/(Ax ) 

chromosphere g ll 


corona 


(32) 
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photo photo 

Wto- 

photo 


(33) 
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Im (Axp j|| . 




(34) 


(35) 


(b) VZAcu&AZon 0^ un ECzctAodijnamic ClncwU. Analogue 

It's clear from equation (28) that the global electrodynamics of the 
magnetic loop system can be modelled via a simple LRC circuit analogue. The 
equivalent circuit elements described by equations (29) *(33) obviously represent 
averages over the loop and its underlying mechanicai driver. Specifically, 

A|| and v^ found in equations (29) and (30) represent a characteristic length 
of the system and a characteristic average value for the Alfven speed vhose 
explicit dependence upon spatial details of the system is at this stage not 
overly Important. In fact, since the majority of the system resides within the 
corona where the Alfven speed is at its maximum value, the expressions for the 
inductance and capacitance can for all practical purposes be written as. 


L 


G 




cor 




(36) 


(37) 


where it is understood that v^ refers to the coronal Alfven speed. An important 
point to be made here is that the magnetic loop naturally supports global structure 
oscillations at a characteristic frequency, 0^, given by, 




o 


1 


(38) 



^/ L C ^cor 

One such example of global structure oscillations which are characterized by 

high shear are Alfvenlc "surface" waves (c,f,» lonson, 1977, 1978; Wentzel, 

1979, and references within these articles). It Is important to 

note that Alfvenlc "surface" waves ^ require a discontinuous Alfven speed 

profile. In fact, any degree of nonunifomlty will support their 

existence. Hollweg (1981) has also stressed the Importance of global loop 

resonances. Thus it appears that the formalism presented In this article casts 

the concept of Qtobal structure oscillations and resonances into a unified 

picture. That Is, according to the derived electrodynamic circuit analogue, 

these global oscillations can be interpreted as glpbal LC oscillations. Of course, 

this point was missed by previous proponents of electrodynamic circuit analogues 

since they, in their purely phenomenological treatment, did not account for a 

critically important energy storage element — the capacitance. The capacitance 

is important because it describes the hydrodynamic reactance of the plasma to 

imposed emfs, an essential feature of all hydromagnetic oscillations. This can 

be simply understood by noting that the energy stored by the capacitor is 
2 

1/2 C I where C * (the free space capacitance across field lines being 

C * i’li/4Tr and K being the plasma dielectric) and Where $ * j 6 E is the average 

electric potential across the field lines. Since the low-frequency dielectric 
~ 2 2 

K * 1 + c /v^, one can see Inmediately that the energy stored by the capacitor 

2 

comprises two parts, viz., the electric field energy per unit volume, E /8TT and 

2 2 2 

the polarization energy per unit volume, c E^/8 ttv^. Noting that 
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the polarization energy li therefore In the fora of kinetic energy* 

l/2pv^» In fact* the so-called polarization current ia defined such 

that X B /c « pv resulting in ■ c^E /4Tlv^, tdiich is consistent 

"*x "*o "*i ® A ! 

pol • ^ 

with a derivation of using particle orbit theory* 

Since the hydrodynamic reactance of the 6 1 plasma (i*c.» related to C) 

is associated with electrodynasiically-induced velocity fields* one must also 

self-consistently include viscous dissipation* Therefore, Just as Joule 

dissipation is an Importance source of resistance* current flow* 

so is viscous dissipation since it tends to prevent the electrodynaoiic charging 

of the 6 < I plasma capacitor by inhibiting the development <nf polarization 

velocity fields within the magnetic loop* Although viscous drag is also important 

within the 6 I photospheric driver* it does not explicitly appear in the 

electrodynamic description of the system. Rather* photospheric viscosity is 

implicitly Included in the driving emf*^(t), which depends upon the photospheric 

velocity field, which in turn depends upon a balance of inertial, non-elec trodynamic 

and viscous forces. Within the 3 < 1 magnetic loop, however, viscous dissipation 

makes an explicit showing in an LRC analogue as a phenomenological viscous 

resistor, R , * Pr R , where the magnetic Prandtl number is defined by eqn. (22) 

viscous mag Joule 

2 

and where the Joule resistance, • Therefore, the 3 < 1 

resistance, given by equation (32) is the sum of Rj^y^g and 

over both the chromosphere and corona* Note that the plasma parameters found 
within the summation of equation (32) represent local conditions with and 
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Ax^ being the field-aligned extent of the region > anAthe local croas-field ahear 
length of the Induced currents and velocity fields which is given by equation (26) 
if local resonance absorption occurs. 

Since the cross-field shear length, represents the cross-field spatial 
extent of the field-aligned current, it is clear that the total field-aligned 
current, Ijj - (Ax^)Jjj, It is important to note that Ij, — a consequence 
of the shear constraint, (7 x 0. Note that cannot possible be zero since it Is 
driven by perpendicular^ polarization E^-fields which are essential for the 
establishment of a non-zero Poynting flux into the loop from its underlying 
driver. Therefore, even though Jj| » seems to indicate that most of the 
induced current is force-free it must be kept in mind that Jjj is confined to flow 

a filament of thickness Ax whereas j flows through a filament of thickness 

r r 

fj^»Ax^. Therefore the total currents, . TT£^(Ax^)Jj| and Ijj - 

are equal as is pointed out in equation (34), It is very important to note 
that steady heating of magnetic loops from an underlying driver is always 
associated with Lorentz forces within the loop which Induce a ^ < 1 velocity 
field. This effect is a natural consequence of the capacitive term in equation 
(28). 

A simple way that the Ax^ can be estimated, which is equivalent to the 
derivation in section Il.a., is to think of the global structure oscillation 
as an exciter of internal field line oscillations within a spatial bandwidth Ax^. 

In other words, there is a mode coupling between the global structure oscillation 



the local eigenfrequency. Noting a correspond encit between the width of a local 
frequency resonance, Aw, and the width of the local spatial resonance, Ax^, i.e., 


Ax 

r 


Ate 

te 


a 


( 39 ) 
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where the Alfven speed *■ gradient scale length » "a" is defined by equation (25). 
Since the width of a frequency resonance depends upon dissipation through an 
oscillators "Q" value» l.c.» 


(i) 



m 


where Q * L/c/R ^the ratio of the oscillator's natural Impedence to the 
resistive impedence, so qiust the width of a spatial resonance depend upon the 
magnetic loop system's local value, i.e,, 


where 



1 

^local 


‘local 


local 


IJ¥) 


local 


(41) 


(42) 


Noting that corresponds to the local resistance described by equations 

(32) and (33), and that the local natural Impedence is proportional to the 
local Alfven speed, equation (42) reduces to, 


Q 


local 



* f'.ag 



(43) 


where the magnetic Reynold's number, Is defined by equation (27) and 

where H(1 > 8) is the Heavyslde step function which is used to Illustrate that 
the viscous resistor does not appear In 8 ^ 1 regions. Using equation (43) in 
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•quAtlon (41) and aolvlng for 4x^ lanodlatoly yloldi oquatlon (26). Thoroforc, 
for a f pacific driving frequency, <e,g,, a loop'e global reeonance frequency, 
0 ^) abiorption will occur along resonant field linea of width given by 
equation (26), provided the condition for the existence of spatial resonances 
is satisfied, i«e«, v^ (x^ ■ x^^^)d/dx|j ^ W, In this regard, the local 
resistances, nomalized to the natural impedence of the entire magnetic loop 
system, l,e,, ^ * *^local * ***P*”‘^ upon the magnetic Reynold's number 

in two different ways depending upon whether or not local resonance absorption 
is occurring, viz., 


local 

R 

local 


+ W.uig “O - 


locals 

^ \ n -1 

^ I “-8 


locally 

non-resonant 


(44) 


(-h) 


-2/3 


i«oiant 

m&g nmg resonAnt 


These normallzud resistances are also related to the global quality, Q . . , , 

global 

of the entire loop system defined as the, ratio of the system's natural 
impedence, V L/C, to the total resistance in the circuit, i.e,. 


^global 


E 

photosphere 

chromosphere 

.corona 


- -1 


local 


(46) 


As I have already mentioned, local resonance absorption does not necessarily 
occur everywhere along the magnetic loop system. However, this process does 
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occur at least in the coronal )>ortion of the magnetic loop where the local 
resonance condition is easily satisfied* 

It is becoming quite apparent that given the driving em£»£(t), which is 
explicitly correlated with the non*elactrodynaaically driven* P ^ I photoapheric 
velocity field (e*f. » equation (35))* it is a simple matter to determine the 
electrodynamic response of the magnetic loop. This respomie comprises both a 
resistive and reactive component. The resistive response is* of course* related 
to the electrodynamic heating function* E * which will be derived in the next 
section. The reactive response is both microscopic as well as macroscopic. 

The .nicroscoplc response is associated wif^h the field-aligned current* J|{, which 
is carried primarily by the electron fluid and which* in the context of the 
one-fluid model used in this article, can be thou^t of as a microscopic ingredient 
of the plasma. The macroscopic response is associated with the cross-field 
polarization current, which is carried primarily by the ion fluid, Tue 
polarisation current* j^* results in a Lorentz force* ^ x B^/c, which drives a 
macroscopic velocity field within the magnetic loop. These reactances will be 
derived in section II. (d) and will yield very important information regarding 
the possible development of microscopic nonlinearities such as anomalous transport 
and macroscopic nonlinearities such as the conversion of sheared flow fields into 
compressible flow fields. 


(c) The GeneAoJUzdd EtixtAcdynmic. HavUng FmcUon, 

In general, the driving emf,^(t), in the global electrodynaalci equation 
(28) ia beet represented as a stationary random variable* This follows from 
the fact that the 8 > 1 velocity field, ^(t), idiich generates the enf,^(t), 
is also a stationary random variable. Kolmogoroff proved in 1940 that every 
stationary random process can be represented by a linear combination of harmonic 
oscillations* If the duration of the process is very long and we 

want good accuracy, one must use an increasingly larger number of harmonic 
oscillations spaced arbitrarily closely. Formally, this is expressed by the 
Fourier-Stieltjes transform. Specifically, the transform of^(t) and its 
inverse are given by 



We can readily appreciate that d^ are like L'ourier coefficients but In this 
representation discontinuous jumps in are allowed. This formalism therefore 
permits analysis of random variables which contain discrete frequencies 
superimposed upon a continuous spectrum, a rremon situation In turbulent 
systems such as the solar convection zone. Furthermore, since we are dealing 
with stationary random processes, it follows that the ensemble average of£(t) and 
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<8(t)>-<dej^>*0 


(49) 


where the breckete repretent en enteoble average. Since the cntcinble average 
of products of random variablee are in general non>zerO| the^.^ are phytically 
tignlficant. Specifically, the spectral "power is represented by 


we=)„«(<ie*de„)=<e5)^du 


(50) 


and is related to <6*(t)S(0)via 


(S*lOe(t)>* J" <d8‘>y 


(51) 


where the * refers to the complex conjugate. 

Armed with this powerful formalism, it is a simple matter to determine 

Che generalized electrodynamic heating function, E which is defined by, 

n 



4 


X cor 


< I*(t)I(t)Rj^_,p > 


(52) 


2 

where TTj£ I /4 is the volume of the coronal loop and wtiere R, is defined 
X cor loop 

by equation (32) which is evaluated by using either equation (44) or (45), 
Taking the Fourier-Stieltjes transform of equation (28) yields. 


6« 


w R -1 u 1 [i + i(<d^w - w /o^>Q - . , ] 
global ' o o ^global 


( 53 ) 
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where ^global defined by equetloni (31) end (46). Noting that 




(54) 


and using this expression along with equation (52) results In a generalized 
electrodynamic heating function which Is given by, 


E 


H 


■V/. 

"glob.1 


where the "power factor." P£(<4 is given by, 


(55) 


Pf(W) 


[1 + («/« - Ui /«) 

0 0 


2 ^2 


^global ^ 


(56) 


The well known plot of p^(a!) versus is illustrated in Figure 2, Clearly, 

for » 1 the power factor is sharply peaked about co « We shall see 

that magnetic loops typically have Q , . , » 1 allowing a simple evaluation of 

global 

2 

the integral in equation (55). Specifically, can be taken outside the 

integral provided it is evaluated at the resonance frequency tti « Ci^. It is then 
a simple matter to evaluate the integral of the power function over frequency, 
resulting in a resonant electrodynamic heating given by. 



( 57 ) 
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vhere 0) • 2lTv has been used. Eliminating Q , . , ■ V L/cXr , . , from 
00 'global / global 

equation (57) yields a result that is virtually independent of dissipation, i.e.. 


' • (^t:) (^) (w) 


(58) 


In fact, the power drain on the driver, E..R , . ,/R, , is apparently totally 

H giODal loop 

independent of dissipation. TfUi h> cl cJiaMJ^lc. izcutufio. oi KUemnee. phenomena and 

doti not -mptij tixat tnAeveutbttitic^ oAe WfUinpoAteuxt, Irreversibilities determine 

2 

the fraction 4v /v ■ TT/o , . , , of the available driving emf, < C >o v , that 
o 0 global ^ ’o o 

interacts with the electrodynamic system at resonance. Since at resonance the 

inductive and capacitive impedences cancel one another, the total power drain on 

2 

the driver equals ^'’o^^global explicitly independent of 

the dissipation since course, the power absorbed by the loop Is simply 

*^loop^^global power drain on the driver (i.e., c.f., equation (58)). 

An important point to note regarding equation (58) is that resonant heating of 
magnetic loops is characterized by an explicit correlation between the spectral 
details of the 3^1 velocity field and the power absorbed by the loop which, 
of course, critically participates in determining the loop's thermodynamic state. 

Noting that 


^2 


. ^ _ - photo 

16TT (vP^°*^° £ /c)^ <1/2 pv-^>v 

A X y ''o 


(59) 


which includes a factor of two stemming from the fact that there are two sources 
of emf, one at each footpoint, equation (58) can be rewritten as 
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H 


16 TT 


( photov y photo ^ ^ 

H (tr-) (&)■ 


2 fhoto 

"''e =’v„''o 


( 60 ) 


In Che context of the above discussion, this result can easily be derived from 
ideal Dihd, noting that all electrodynanic energy entering the loop will be in 
reality dissipated — a point stressed by lonson (1977, 1978) who argues in favor 
of an underdamped, high-quality resonant system with also 

stressed by Rosner et (1978) and Golub et «1 . (1980) who argue in favor of 
an overdaraped, low quality system with Specifically, the net flux 

of energy into a loop is the product of the transmittance, 2v^ input 

signal speed, the energy density of the 8^1 photospheric velocity field, 

< 1/2 pVg number of polarization states which equals two for 

electrodynamic disturbances, the number of normal modes wliich equals two for 
low-frequency electrodynamic disturbances in 8 < 1 plasmas, the number of sources 
which equals two for a magnetic loop with two footpoints, and the fractional 

dissipation within the loop, Rioop/*^global* power absorbed 

per unit volume is the divergence of the above product of terms with V -* 
one readily recovers the generalized electrodynamic heating function, E^, given 
by equation (60). Since « >J 3hgT^^^/m^ where is the 

effective black-body temperature of the star, the heating function can also be 
written as , 


E 


H 


6.28 X 10 



loon 
^global 


photo 

< 1/2 pv, 



( 61 ) 
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Equations (60) or (61) are of major importance since they clearly reveal 
how the electrodynamic heating function scales with the spectral power of the 
6^1 photospheric velocity field (i*e., goal (i) of the Introduction). The. 

&6i&ntiixZ (JcuLtuAc oo a Ae.iona'it heating mecfuintdmi^ that magnetia toop& uUth cUi^eaent 

Zenciths/'and hence di^^eAent gZobat aeionance iAequenciei, 

one heated c<t a Kate that cAiticattij depends upon the amount o^ ^p^etKot 

poioex ^ound at vtz. , < 1/2 pv^ This is in complete contrast to a 

non-resonant heating mechanism (i.e., Q . . . < 1) which would depend upon the 

global 

total power (i.e., the spectral power summed over all frequencies) rather than 

the details of the photospheric velocity field's spectral power function. Thus, 

although all electrodynamic energy entering the magnetic loop would be iinnedlately 

dissipated by a non-resonant heating mechanism, as has been pointed out by 

Rosner et ajL. (1978) and Golub et al. (1980) there would be no correlation between 

the amount of energy absorbed by a loop of length (and hence its temperature 

and pressure) and the photospheric spectral power at » '^A^^'^cor* also 

that although the role of the magnetic field, B , appears to be passive for a 

o 

resonant heating mechanist^ in the sense that it part'clpates in determining a 

loop's global resonance frequency, v * v /2i , the photospheric power spectrum, 

o A cor 

< 1/2 pVg could in itself actively depend upon B^. In this regard, 

independent investigations of the role of magnetic fields in 6 ^ 1 turbulence will 
be invaluable in shedding light upon the nature of the photospheric velocity 
field's spectral power as it appears in equations (60) and (61). 
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((f) The. ReacUvz Ete^tAodynamic heipontz and the. Hole ojJ jAAzvzulbltttieA 

since the nejority of the megnetlc loop eystem Is found within the corona, 

this section will focus upon the reactive response of the coAonaZ portion of nagnetlc 
loops to electrodynanlc driving. Furthermore, I will assume that a local resonance 


absorption process is occurring within the corona since the local resonance 
condition is easily satisfied. 

As I have already discussed, the power drain on the photospherlc driver is 
explicitly independent of the dlssipatiovi -- a natural feature of resonance 
phenomena. This, however, does not mean that irreversibilities completely 
disappear from the problem. In fact, the root>mean square amplitude of the 
induced current, idiich is associated with the microscopic (l.e. , electronic) and 
macroscopic (l.e., ionic) reactance of the plasma contained by the magnetic loop, 
does indeed depend upon irreversible processes such as electron>ion collisions. 


The microscopic reactance is related to the field-aligned current density, , 
which is carried by the electrons drifting at velocity v^ along the ambient field. 

The macroscopic reactance is related to the cross-field current density, , which 
through Lorentz forces drives an azimuthal flow of coronal plasma with velocity Vg. 
Noting that the dominant drag on the ions is viscous drag whereas for electrons it is 
Joule drag, it follows that the ratio of ion to electron heating is the ratio of the 
viscous resistance to the Joule resistance, l.e.. 


viscous ion heating 
Joule electron heating 


viscous 


R 


Joule 


Pr 

mag 



( 62 ) 


where represents viscous drag on the azimuthal flow of plasma and 

where represents ’’Joule" drag on the field^allgned flow of 


electrons (c.£., section II. (b)). Again, there appears to be no explicit 

dependence upon Irreversibilities such as collision frequencies. The appearance 

of Irreversible phenomena Is, however, explicit In expressions for the electron 

drift speed, v^, the plasma's azimuthal flow velocity, v^, and the width 

of the coronal field lines which are subject to electrodynamlc activity. 

The width (Ax be easily estimated by noting that a magnetic loop 

* 

responds at frequencies other than Its global resonance frequency, <i^, given by 
equation (38). In fact the loop will respond for a continuum of different 

frequencies, Acij^ ~ '^o^^global* *^global global quality of the entire 

loop system defined by equation (46). As one might expect, Aw^ results In a 
continuum of excited spatial resonances whose overall thickness (4x^)^^^> is given by 


Ax 


(Ax ) » •— 

r cor R 


cor ^global 


(63) 


where it is understood that ^ ^ represents the local quality of the coronal 

cor 

portion of the magnetic loop given by equation (45) and where In this case, Ax^ 
represents the width of a single spatial resonance within the corona given by 
equation (26). 

The electron drift speed, v.. Is estimated by noting that, 

d 


1 
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R 


it Vc/l 


global 


(64) 


29 


with I • ■ Ijj itaa equation (34). Since Jjj ■ nev^, it limply 

follows from equation (64) that the electron drift speed normelieed to the ion 
sound speed is given by. 




(65) 


photo 


where equations (26) and (59) have been used. In equation (65) v^ represents 

the ion Sound speed, B . ^ /B and T .JT represents the ratio of the photospheric 

pnoto o ex£ 

magnetic field to the coronal field and the ratio of the photospheric temperature 
to the coronal loop temperature, represents the ratio of the plasma skin 


depth (i.e. , Alfven speed's 


scale length (c.f. , 


equation (25)), ^ is the normalized coronal resistance given by equation (45), 
cor 

Q , . . is the global quality of the magnetic loop system defined by equation (46) 
global 

and (Vb) . ^ represents the velocity of the photospheric driver (i.e. 
o pnot<) 

Vq ~ J ' 


assume that the coronal resistance dominates both 


the chromospheric and photospheric resistance -- a situation which occurs for 

solar applications (c.f., section III), then Q , . , ^ and equation (65) 

global cor 

can be written as 


(■ 
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( 66 ) 


photo 


where equation (45) for R has been used. For the case of solar magnetic 

cor 

loops, equation (66) can be further simplified by noting that “ 3.38, 


I 


30 




solar 


cor 


4.2 X 10‘® (I + Pr 
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photo 


(67) 


The azimuthal flow velocity, i> estimated by noting that, 


1/2 C W “• Q , . , I^R , . , 
o ^global global 


( 68 ) 


where $ (v.) B 4 /c is the coronal loop’s capacitive emf. Equation (68) 

t7 cor o X 

results In the following expression for the coronal flow velocity, normalized to 
the coronal ion sound speed; 


i-i-lr ■ i-^) i-H 


1/2 


1/2 

^global 
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(69) 


photo 


which, for the solar parameters used earlier, reduces to 


(■?■) 


solar 


cor 


0.25 (1 + Pr )'^^^ Re^^® 

mag mag 


(t) 


(70) 


photo 


The size of the solar coronal region within which the resistive and reactive 

activity described by equations (61), (66) and (70) occurs is given by equation 

(63) for (4x ) 

' r cor 


(Ax ) 

r Gor 

jt 


(l^Pr 


mag' 


mag 


( 71 ) 


where equation (26) for the width of a single resonance, 4x^, has been used 

along with Q , , . 

global cor 


(c) NonlincaA. E^iectd 


Knowledge of the plaima reactance derived by a linear fornallim la very 

useful In determining the relevance of various nonlinear processes. For 

example, microscopic nonllnearltiss associated with the destabllzatlon of 

the field-aligned electron current could result In an anomalous Increase In 

the electron collision frequency, » Slven by equation (14). Such 

microscopic nonlinearities can only occur if (v./v ) >1 and would result 

a 8 cor 

in both anomalous current dissipation and anomalous viscous dissipations 

Since the magnetic Prandtl number, Pr , as derived from viscous dissipation 

mag 

of sheared, incompressible flows^ depends upon the Ion collision frequency wliich 

in turn depends upon the electron collision frequency (c.f. , equation 

anomalous 


(11)), it follows that the anomalous Prandtl number, Pr 


mag 


, is given by 


Pr 


mag 


Pr 


anor:talous 

mag 


(•^r (‘ • (^) (t)f 


(72) 


where v , has been used (i.e., equation (14)). Thus the ratio of 

anomalous pi 

ion to electron heating described by equation (62) could be significantly 

affected by the onset of microscopic plasma turbulence. In addition, the 

themal vs. nonthermal energy partitioning of the "heated” electrons is affected, 

with the nonthermal component becoming increasingly enhanced as the ratio 

(v,/v ) increases beyond unity (c.f., lonson, 1981). The width of the 
d s cor 

spatial resonance, 4x^, would also be increased since the magnetic Reynold's 
number. Re , depends upon the electron collision frequency. In fact, 
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would decrease with the development of microscopic plasma turbulence also resulting 
In a decrease in the macroscopic flow velocity* (v») » described by equation 

(70). 

A second possible nonlinearity that could develop would be driven by the 
macroscopic flow field, Although for a linear analysis this Is a sheared flow, i 

compressible component could be generated if I" this case it 

would be necessary to consider viscous heating of compressible flows which is 

much more efficient than viscous heating of sheared flows (c.f., Hollweg, 1979). 

In fact the relevant coefficient of viscosity used in deriving the magnetic 

Prandtl number would be Hy rather than u . Therefore, if (v./v ) ~ 1 then 

11 X ' 0 s cor 

the magnetic Prandtl number, py®^P*^®*®ihle ~ g(nj^/in^)^^^ (0^/v^). In general 

2 

however, the flow is not completely compressible. Since (v./v ) is a measure 

0 8 cor 

of compressibility, the effects of compressibility can be approximated by using 
Che following expression for the Prandtl number: - 

This could also significantly affect the ratio of ion to electron heating as well 

as Che width of the spatial resonance, Ax . 

r 

It is quite clear that the ion to electron heating rate depends critically 
upon the plasma conditions of the system of interest as well as whether or not 

I 

microscopic and/or macroscopic nonlinearities develop. Note that if both types I 
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of nonlinearitlet occur, then the anoMlout Ion colliiion frequency, 

1/2 

should be used in equation (73). It is iaportant to 
note that the question of nonlinearities can only be properly addressed by 
first investigating the linear response -> specifically, the reactance 

discussed in section II. d. It will become apparent in the next section that applications 

to solar coronal loops will require nonlinear modifications. 

III. Solar Applications 

In order to apply the physics of electrodynamic coupling described in 
section II., it is necessary that we have some knowledge of the magnetothermo- 
dynamic conditions characteristic of the particular plasma setting of interest. 

As a cimon-itniitLon of how this physics can be utilized, a specific application 
to solar magnetic loop systems will be presented. 

(a) Ro,ionaivt dynamic. Heating oi SoZoa CofionaJL Loop6 

Serio al. (1980) have recently generalized the scaling laws derived 
by Rosner et al. (1978) to include solar coronal loops whose height exceeds 
their pressure scale height. These scaling laws, which are fairly consistent 
with observations are given by 


T 


1.4 X 10“ 


cor 


1/3 


exp 


E 


loss 


,„5 „1.17 

10 p ' 


,-0.83 

'cor 



(74) 


(75) 


where T is the maximum temperature of the coronal loop, P the base pressure, 
3 

6.12 X 10 T cm, the pressure scale height and where 


is the average 
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energy per unit volume thet is lost through chromospheric end coronsl rsdlstlon. 
In equfitions (74) end (75) it has been assumed that the energy deposition scale 
height is large compared to the loop length since we are dealing with a global 

a 

heating mechanism* 

Golub et al. (1980) have also derived an empirical relationship between 

the base pressure of coronal loops and the average strength of the underlying 

photospheric magnetic field* Assuming that the coronal magnetic field* 

B * can be estimated by a simple dipole extrapolation of B . . (i*e.« 

o pnoco 

3 

B /B . ^ (2<^ f2l ■¥ t ) “" 0.296 with a dipole displacement of 

o photo coir cor cor 

2X^^^)* their relationship becomes* 
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o , 1.6 

2.1 X 10 B 

o 



(76) 


From equations (74) and (76) it is possible to determine the magnetothemodynamic 

state of a coronal loop as a function of its maximum temperature, T* and length, 

a . Specifically* 
cor 


p 


3.64 X 10 


-10 


B * 1.41 X lO’^ 

Q 



(77) 


(78) 


and, noting an exponential drop in the coronal density with height, the coronal 


Alfven speed is given by. 
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2.68 X 10 



(79) 


mmn * 

In order that the tefflperature» Ti^be eliminated ^ Equations (77) • (79) It la 
necesaary to utilize equation (75) for noting that 


E ■ E. 

H loss 


(80) 


where E„ Is the generalized electrodynamlc heating function given by equation 

(55). A simultaneous solution of equations (77) - (80) will then result in a 

scaling of loop temperature. T. and base pressure. P. with the loop size, 

4 , and the spectral power of the photospherlc velocity field, 

cor* X 

< 1/2 evg (l.e», goal (11) as stated In the Introduction). 

To determine the relevant form of electrodynamlc heating (l.e. . resonant 
or nonresonant) it is first necessary to estimate the global quality. 
of the loops. This can be easily done by using the following parameters as 
&&tancutQA> in equations (27). (44) - (46). 
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The rcvults ere: 


l^ePhoto ^ ^ ^ ^ ^ 

mee photo 


-9 


Ke"^'“ ~4»l0’i , -2x10-’ 

mag mag "‘chrom 


„ cor - ,-12 „ cor , ~ „ ,, 

Re ~ 2 X 10 ; Pr ~ 4.3; R 2 x 10 

mag mag cor 


^ L/C 0,24 Ohms 


Q , , , ~r”^ 5 X 10^ 

^global cor 


where it has bee» assumed that local resonance absorption occurs within the 
corona (l.e., equation (45) for was used to estimate . The rationale 
for assuming tha*; local resonance absorption occurs within the corona is based 
upon the fact that the global resonance frequency of the entire magnetic loop 
system depends primarily upon coronal rather than chromospheric and photos pheric 
conditions (c.f,, equation 38)). 

It is quite obvious that the global quality of typical solar magnetic 

3 

loops is extremely large (i.e., ~ 3 x 10 ) -- at least under the 

constraints of a linear theory. Therefore, the resonant form of the heating 


function, E^, described by equation (61) should be used in the energy balance 
described by equations (75) and (80). This results In the following /..mpllclt 
equations for the aMxlsiuin temperature, T, and base pressure, P: 


3.36 X 10^ [t ,, exp (0.36je H ) < l/2pv? 

eff ■ cor p ^0 


photo 1 
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. exp / 1.51 I ff < pvg 

X \ cor p / eff 9 

cor \ ^ J _ o 


, 6/7 


^dynes/cm^^ (82) 


idiere R, R . i. , has been used, 

loop global 

In order to solve equations (81) and (82) explicitly for 2 and P as a 

function of loop length, It^ Is necessary that we know the spectral 

power of the photospherlc velocity field. As a demonstration, the power 

spectrum Illustrated In Figure 3e has been assumed. This spectrum corresponds 

to observations of solar plage regions where densities are of the order of 
-7 3 

10 gm/cm . Illustrated are two major peaks, one at Intermediate periods 
(300 > r > 100 sec) (Orrall, 1966; Woods and Cram, 1981) and one at longer 
periods (r > 400 sec) (Harvey, 1980). Figure 3d Illustrates how the global 
resonance period versus loop length scaling, the spectral 

power function of the photospherlc velocity field into 

equation (61) (c.f.. Figure 3c)^T » equation (81) (c.f.. Figure 3b) 

and P “ equation (82). Also plotted in Figures 3a and 3b are 

observations of active region loops (Landinl et_ al . , 1975; Pye et al . , 1977) 
wlilch absorb power from the "Intermediate regime" of the photospherlc spectral 
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power function; large-scale loop structures (Maxson and Valana, 1977) which 

are electrodynamically coupled to the "local power minimum" region of the 

photospherlc spectral power function; and very large scale active loops (Pye al, , 1977) 

which can absorb power from the "long period" regime of the photospherlc spectral 

power function. 

(fa) The Reactive Res pome and the Rote 0(( \HAevQM.ihiZitits tn SotoA OoKonaZ Loops 

Using the results of Fig. 3. Fig. A illustrate^ how the natural Impedence, ^ L/C^, 

magnetic Prandtl number, Pr , and magnetic Reynold's number. Re of solar 

mag mag 

. -2 

magnetic loops scales with the loop length, £ . Noting that (Vq/v ) . 5 x 10 , 

cor o 8 photo 

these parameters are used In equations (45), (46), (67), (70) and (71) for 

Q 1 u i» (^j/v ) » (Va/v ) and(Ax ) which were derived from a linear analysis, i 

^global’ d s cor 6 s cor r cor •' * 
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which yield the results illustrsted in Figure 5. One noteworthy feature of 

this illustration is the extremely large Q-values leading one to believe that 

the microscopic and/or macroscopic reactance stenming from a linear theory could 

be quite appreciable. It is Important to note from this figure that although 

microscopic nonlinearities associated with anomalous current dissipation are 

LUvimpoAtant (i.e., because (v./v ) « 1) macroscopic nonlinearities associated 

a 8 cor 

with the conversion of a sheared^ incompressible coronal flow into a partially 

compressible flow probably are important because ~ This is 

consistent with Hollweg's (1981) independent analysis of the problem. Therefore, 

compressible 

as discussed in section II. e., equation (73) for Pr should be used. 

mag 

This will result in increased viscous dissipation which lowers the Q-value of 
the loops^ thereby resulting in both a broadened spatial resonance and a decrease 
in the coronal flow velocity, 


nonlinear _ , ^ “1/12 1/4 

Q , ^ , =“ 1.6 X 10 Pr Re 

^global mag mag 


(87) 




nonlinear 


-3 - 1/8 1/8 
2.2 X 10 Pr Re ' 

mag mag 


( 88 ) 


cor 


,, ^nonlinear . .nonlinear 

(Ax ) =“ 0.1 X /Q , , - 

r cor global 

cor 


(89) 


Figure 6 Illustrates the nonlinear reactance and Q-values described by equations 

(87) • (89) for the case of solar magnetic loops. Note that although the quality, 

nonlinear^ is significantly smaller than that estimated from a purely linear 
global 

analysis, it is still larger than unity validating 


use of the resonant 
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heating function, E . Thli nutt be to alnce It £• the reeonent charecterlstics 

M 

of the system uhich nonllneecly supports the compress ibllity of the flow. 

Therefore, the results depicted in Figures 3 and 4 remain unchanged. 

As is evident in Figure 6a, the cross -field size of the coronal dissipation 
site has dramatically increased from the order of kilometers (c.f.. Figure Sa 
and lonson, 1977, 1978) to the order of thoaiandi of kilometers due to the 
nonlinear modifications presented here. This is an extremely important result 
since there has been a great deal of concern regarding how solar coronal loops 
can be heated throughout their entire volume (c.f., Hollweg, 1979). 

Another interesting feature of this analysis is the prediction of an 11-16 km/sec 
coronal velocity field (c.f.. Figure 6b) and Che preferential viscous heating of Che 
ions, l.e., 

viscous ion heating ^ . compressible ^ 8 

Joule electrcr. heating mag ' ' 

Both of these effects could be responsfble for the widespread observations of 
coronal line broadening (c.f., Feldman and Behring, 1974; Acton et al», 1981). 


IV. Conclusion 

The major emphasis of this article has been the physics of global electrodynamic 

coupling between a 0<1 magnetic loop and an underlying 0>1 mechanical energy reservoir, 

A rigorous analysis of this problem has revealed that the physics can be represented 

by a simple yet equivalent LRC circuit analogue. This analogue points to the 

existence of global structure oscillations (i.e., v «v . /2A ) which resonantly 

o A cor ^ 





PLASMA FLOW THICKNESS OF NONUNEAR 

VELOCITY. SPATIAL RESONANCE 


PIG 6 



excite Internel field line oeclllatlons at a apatlal resonance within the magnetic 


loop (l.e., Ax^ being the width of the spatial resonance). Although the width of 

this spatial resonance as well as the Induced currents, and explicitly depend 

upon viscosity and resistivity, the ACAonant form of the generalized electrodynamic 

heating function, E^, Is virtually Independent of Irreversibilities. This Is a 

classic feature of high quality resonators that are externally driven by a broad 

band source of spectral power. This follows from the fact that at resonance the 

heating function depends solely upon the resistance, R, and the emf which can 

intZAAoX. with the loop, <r >,Av , l.e., £„*< C >,,Av /R. Since the interaction 

^ o li V o 

° ° r I 

bandwidth, Av^, also depends upon the resistance, l.e*, Av^/v^^l/Q • rJc/L, it 

follows that Ejj“< £ > V^/J^ L/C* which is explicitly independent of R, In addition, 

2 ^ 

since < t “<^pv It becomes clear how the heating function scales with the 
o o 

spectral power function of the^ mechanical driver. The essential feature of a resonant 
heating mechanism is that magnetic loops with different lengths and hence different 
global resonance frequencies are heated at a rate that critically depends upon the 


amount of 3^1 spectral power at the resonance frequency, V^. 

As a demonstration, this physics was applied to the solar setting. A linear 
analysis resulted in extremely large Q-values for solar coronal loops Implying very 
narrow spatial resonances of the order of kilometers and a coronal flow velocity 
of the order of the sound speed. The electron drift velocity was found to be orders 
of magnitude less than the ion sound speed and therefore microscopic plasma instabilities 
were not important. However, in light of the large coronal flow velocity it was 
necessary to consider nonlinear modifications to the linear analysis stemming from 
the development of a compressible component to the coronal flow. The nonlinearly 
generated compressible component of the flow allows the use of a viscosity coefficient 
which dramatically Increases the rate of viscous dissipation. As such, the steady 
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state quality of solar coronal loops Is decreased to Q««15 which results In broad 
spatial resonances of the order of thousands of kilometers » a coronal velocity field 
of the order of 10-16 km/sec and preferential viscous heating of the lonst The 
nonlinear modifications, however, do not change the temperature (pressure) scaling 
of solar coronal loops with their length. This follows from the fact that the system 
Is still highly resonant (l.e>, Q>1) allowing one to utilize the resonant form of 
the heating function derived from a linear analysis. It Is remarkable that such a 
simple theory Is so consistent with a variety of observational constraints, despite 
the neglect of nonlinearltles other than those addressed In this article. 
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Prototype ugnetlc loop systen and equivalent LRC circuit* 

Power factor! p^ (<>>)> versus the frequency nomalized to the resoance 
frequency! 

Illustration of how the solar photospheric pov;er spectrum eloctrodynamically 
couples to coronal loops. The data points refer to active region loops 
(l.e., X's ; Land ini et al.. 1975; Pye^al.! 1977>! large-scale structures 
(l.e., crosses; Maxson and Vaiana! 1977) and very large scale active loops 
(i.e., circles; Pye^ aL,! 1977). 

Characteristic properties of solar coronal loops with -0.1. 

cor 

Results of a linear analysis applied to solar coronal loops. 

Results of a nonlinear analysis applied to solar coronal loops. 
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